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ABSTRACT 

A finite-volume lower-upper (LU) implicit scheme 
is used to simulate an inviscid flow in a turbine cas- 
cade. This approximate factorization scheme requires 
only the inversion of sparse lower and upper triangular 
matrices, which can be done efficiently without exten- 
sive storage. As an implicit scheme it allows a large 
time step to reach the steady state. An interactive 
grid generation program (TURBO), which is being deve- 
loped, is used to generate grids. This program uses 
the control point form of algebraic grid generation 
which uses a sparse collection of control points from 
which the shape and position of coordinate curves can 
be adjusted. A distinct advantage of TURBO compared 
with other grid generation programs is that it allows 
the easy change of local mesh structure without 
affecting the grid outside the domain of dependence. 
Sample grids are generated by TURBO for a compressor 
rotor blade and a turbine cascade. The turbine cascade 
flow is simulated by using the LU implicit scheme on 
the grid generated by TURBO. 

NOMENCLATURE 

A , B Jacobian matrices 

Cjj a control point 

0 difference operator 

Ej(r> control point curve 

e total energy/volume 

F flux vector in the x-direction 

Fj(t) control point curve 

f flux vector in ^-direction 

G flux vector in y-direction 

G a integrated interpolation function 


g flux vector in n-direction 

H|} integrated interpolation function 

J Jacobian of coordinate transformation 

M number of control points in ^-direction 

N number of control points in n-direction 

P(r,t) position vector 

p static pressure 

P t total pressure 

Q(r , t) Boolean sum transformation 

R residual 

r curve parameterization 

r A »r B spectral radi i 

T(r , t) tensor product transformation 
T r transpose of matrix 

t curve parametrization 

U,V contravariant velocities 

u , v , Cartesian velocities in x- and y-directions 

W vector of flow variables 

x,y Cartesian coordinates 

on-off boundary switch 

(3 numerical parameter between 0 and 1 

Y ratio of specific heat 


1 



5,n transformed coordinate system 

X eigenvalues of Jacobian matrices 

p density 

x time 

v|/ interpolation function 

INTRODUCTION 

An unconditionally stable implicit scheme that 
has error terms* at most, of the order (Ax) 2 in any 
number of space dimensions can be derived by lower- 
upper < LU) factorization (Jameson and Turkel, 1981). 
Splitting ensures the diagonal dominance of lower and 
upper factors and makes use of the built-in implicit 
dissipation. The LU implicit scheme requires only the 
inversion of sparse triangular matrices* which can be 
done efficiently without using extensive storage. As an 
implicit scheme it allows a large time step to reach 
the steady state. This approximate factorization 
scheme was demonstrated to be robust and efficient in a 
broad flow regime; for example, Jameson and Yoon (1987) 
for transonic flows and Yoon and Jameson (1986) for a 
high speed inlet flow. These previous efforts focused 
on improving efficiency of the LU implicit scheme by 
using multiple grids with simple isolated airfoils and 
inlet geometries. In the work presented herein the LU 
scheme is applied to obtain the solution of the invis- 
cid compressible flow in a turbine cascade. This is a 
first step towards simulation of flows in turbomachinery 
that have more complex phenomena. One such example 
could be a flow in a supersonic through-flow fan 
(Schmidt et al . , 1987). Eventually, for viscous calcu- 
lations the meshes will need to be very fine to resolve 
the boundary layer. And it is likely that the time 
step imposed by an explicit stability bound will be 
much less than that imposed by the accuracy bound of an 
implicit scheme. Since an obvious way to accelerate 
convergence to a steady state is to increase the size 
of the time step, an implicit scheme is expected to 
have a faster convergence. Although the alternating 
direction implicit (ADI) scheme has been valuable in 
two-dimensional problems (8eam and Warming, 1978), its 
inherent limitations in three dimensions suggest an 
alternative approach. 


tions. Eiseman (1987) demonstrated how mesh structure 
can be changed with control points by using simple duct- 
like geometries. The control point formulation is used 
herein to generate grids for turbomachinery geometry. 

A menu-driven interactive program called TURBO Is 
being developed to generate wel 1-structured grids for 
turbomachinery flow field simulations. By using the 
control point form mentioned previously, TURBO can 
achieve a balance between the constraints of the turbo- 
machinery geometry and the desired grid structure. 

A distinct advantage of TURBO compared with other grid- 
generation programs is that it allows the easy change 
of local mesh structure without affecting the grid out- 
side the domain of dependence. The boundary surface 
can be either a fixed form (transf ini te) or a free 
form, which can be changed. Once the grid is improved, 
it can be closely examined by zooming, translating, and 
rotating. Further improvement to the grid structure 
can be made with this closer view of the grid. The pro- 
gram, TURBO, is being built around the control net 
which guides the coordinate curves. By choosing menus, 
TURBO now normalizes the grid near the vane or blade 
surface, provides slope continuity across the periodic 
boundary in C-type grids, stretches grids by stretching 
the control points in real time with the movement of 
the workstation mouse, and simultaneously translates 
multiple control points on a line in a manner similar 
to points on a rubber band. These are some of the fea- 
tures that provide convenience to users. 

The finite-volume LU Implicit scheme was used to 
simulate an inviscid flow in a turbine cascade on a 
C-grid generated by TURBO. The grid generated for this 
flow is normal to the vane surface and has slope conti- 
nuity across the periodic boundary. For this cambered 
blade, the rubberbanding feature was used to reduce the 
grid skewness on the suction side of the vane. Com- 
puted results are compared with the experimental data 
obtained by Goldman and Seasholz (1982). 

For clarity and completeness of this paper, govern- 
ing equations, the LU approximate factorization scheme, 
and the control point formulation are briefly presented 
in the following sections, even though detailed deriva- 
tions are available in the references cited above. 

GOVERNING EQUATIONS 


Flow passages in turbomachinery are highly con- 
strained geometric regions. Fine meshes are needed to 
capture flow phenomena where flow properties change 
rapidly in several subregions. The orthogonality 
desired in highly cambered blades conflicts with the 
periodicity constraints that are required by many flow 
solvers; the need to cluster the grid conflicts with 
the smoothness desired. Generating grids for turboma- 
chinery, therefore, requires a sensible compromise 
between geometric constraints and desired grid struc- 
ture. One of the best approaches would be to generate 
the grid interactively by using a method that allows 
easy local control over the grid distribution (Choo 
et al . , 1988). The control point form of algebraic 
grid generation was formulated by Eiseman (1987). From 
the control points sparsely distributed over the flow 
domain, the shape and position of coordinate curves can 
be adjusted while the grid conforms precisely to all 
boundaries. As an algebraic method, the control point 
form provides explicit control of grid structure and 
requires relatively few computations. The basic ele- 
ments of this method are the multi surface transformation 
for control in a given direction (Eiseman, 1982) and 
the Boolean sum operation for the combination direc- 


The conservative form of the Euler equations in 
Cartesian coordinates for two space dimensions is 


3W 8F 3G 
3x + 3x + 3y 


0 


(1) 


where H is the vector of dependent variables and F 
and G are convective flux vectors as follows: 


W - (p, pu, pv, e) r (2a) 


F = [pu, pu 2 + p, puv, u(e + p)j r (2b) 


G * [pv, puv, pv 2 + p, v(e + p)j r (2c) 

where p, u, v, e, and p are density, Cartesian veloc- 
ity components in the x- and y-di rections , total energy, 
and pressure respectively. The total energy e is the 
sum of the internal energy and the kinetic energy per 
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unit volume of fluid. The system of Eqs. (1) and (2) 
is solved with time x for the primary unknowns, 

(p, pu, pv, and e), and the pressure is obtained subse- 
quently by the equation of state as, 


Here, 0^ and D n are central difference operators 
that approximate 8/95 and 9/8n. An LU factorization 
scheme is introduced to approximate Eq. (7) as 


p = ( y - 1 ) I e 


<u 2 + v 2 ) 
P J 


(3) 


where y is the ratio of specific heats. Because of 
the constrained flow passage geometry of turbomachi nery , 
a coordinate transformation, is performed in such a way 
that the transformed coordinate coincides with the flow 
boundary. With the coordinate transformation, 5 = 

5(x,y) and n = n<x,y), the vector form of the Carte- 
sian components of Euler equations is 


[i + fiJ At(d"A + + D~B + )] 

[i + |3J At(d*A~ + D*B')] 5W + At R = 0 (8) 

where and D~ are backward difference operators 
and Dc and D* are forward difference operators. 
Here, A + , A", B + , and B“ are constructed so that the 
eigenvalues of M +" matrices are nonnegative and those 
of matrices are nonpositive, that is 


9W 

9x 



(4) 


where f * (£ X F + 5yG)/J, g = (nxF + nyG)/J, and 
J = 5 x ny - 5ynx- With the relations £ x = Jy^, 
nx = -Jy^. 5 y = -Jx n , ny = Jx^ and the use of con- 
travariant velocities, 0 = y n u - x n v and 
V * -y^ru + x^v, vectors f and g in Eq. (4) are 
written as 


where 


A+ - i ( A + '•a 1 )’ 

8 + . i (B + r B l). 

r A i max (l X A l)> 



(*- 

r A 

B- .1 

0 - 

r B 


r g > max(|X B |) 


pU 

puU + y n p 
pvU - x n p 
\<e + p)U, 
LOWER-UPPER IMPLICIT SCHEME 


pV 

puV - y^p ' 
pvV + x^p 
V(e + p)V> 


(5) 


To advance the Euler equations with time we 
discretize Eq. (4) in time as 


H n+1 - H n 
At 


KH + IS) + (1 - p)J (i{ + I?) ’ 0 


( 6 ) 


where (3 varies between 0 and 1. For [3 = 0, Eq. ( 6 ) 
is a purely explicit scheme, and for (3 = 1/2, a second 
order accuracy holds in time. Nonlinear terms at time 
level n + 1 are linearized by the following 
approximations: 




A n 6W, 



©" * B " 5H 


where SW = W n+ 1 - W n , and A and B are Jacobian 
matrices defined to be J~^ 3f/3W, and J -1 8g/9W, 
respectively. The truncation error of the above expres- 
sion is of 0 (At 2 ). The linearized implicit scheme for 
Eq. ( 6 ) can be formulated as 


Here \/\ and Xg represent eigenvalues of Jacobian 
matrices. Equation ( 8 ) is inverted in two consecutive 
steps. The LU implicit scheme requires the inversion 
of sparse triangular matrices, which can be done effici- 
ently without using extensive storage. This scheme has 
only two factors in three space dimensions. 

To make the central difference scheme (Eq. ( 8 )) 
stable, an artificial dissipation term in the fourth 
order difference form is added explicitly to R in 
Eq. ( 8 ) for the entire flow region. For transonic or 
supersonic calculation, additional dissipation needs to 
be added implicitly in the neighborhood of shocks 
(Jameson, 1987). In the present subsonic case only the 
fourth order explicit dissipation is employed. 

The cascade flow measured by Goldman and Seasholz 
(1982) is fully subsonic. To simulate this flow, the 
far-field boundary conditions used for an isolated 
airfoil by Jameson and Yoon (1987) were replaced by 
in-and-out flow boundary conditions and periodic 
boundary on the sides. At the inlet boundary, three 
out of four flow variables were specified and the 
remaining one was obtained as a solution of the flow. 

At the exit boundary, only one flow variable was 
specified. In the present work, the density p was 
allowed to change while other variables were fixed at 
the inlet boundary. At the exit, only the static 
pressure was specified. Periodic boundary condition is 
enforced on the sides because of the cascade configura- 
tion of the present problem. 

CONTROL POINT FORMULATION 


[i + |3J At(d^A f D^jsW + Ax R = 0 (7) 

where I is the identity matrix and R is the steady 
state part of Eq. (4) given as 



D c f n + D g n 
5 n* 


The control point array is a sparse grid-type 
arrangement of locations in physical space with an 
index for each direction. In two dimensions it will be 
denoted by (Cjj). As an algebraic method, the control 
point form provides explicit control of the physical 
grid shape and spacing through the dynamic movement of 
the control points. Figure 1 shows an example of a 
two-dimensional control point array. 
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A fundamental part of the control point formulation 
is the construction of curves. This construction repre- 
sents algebraic coordinate generation in a single direc- 
tion wherein two opposing boundaries are connected by 
the newly created curves (e.g., the curve £ 2 ^ 

Fig. 1). With the restriction to only a single curve, 
the opposing boundaries are each represented by a point. 
The first and last points lie on opposing boundaries and 
are the fixed end points of the curve. The remaining 
points are in the interior of the sequence and are used 
to control the shape of the curve. To enforce the 
successive assumptions of the desired direction in a 
smooth manner, a continuous direction field is obtained 
by interpolation. The independent variable for the 
interpolation is simply the curve parameterization. 
Altogether, the interpolated result defines the field 
of vectors that are tangent to the desired curve and 
is simply stated as an interpolation of the first para- 
metric derivatives. Thus a smooth first derivative of 
the entire curve is determined. The desired curve is 
then obtained by a parametric integration. 

Let Cij, C 2 j. .... C^j. be the given sequence of 
M points inspace (e.g., C] 2 , C 22 > C 52 for 5*2 

in Fig. 1); let r be the curve parameterization; let 
Ej(r) be the position at r along the desired curve; 
let n, r£, . .., r M _] be the successive parametric 
locations to interpolate the directions of (C 21 - C]j), 

<C 3 j - C 2 1 > <C M j - C M _1 3); and let Y,, Y 2 

*fM-l be the corresponding interpolation functions which 
successively separate each direction by assuming a non- 
zero value at the associated location while vanishing 
at the remaining locations for interpolation. With 
this notation the desired curve is given by 

N-l 

E j (r) * C U + 2 G a <r) [ C a + l,j ' C aj] 

a= 1 

j*l,2, M (9) 


where 


G (r) 
a 


t a <p)du 


(10) 


To apply the method, the interpolation* functions 
must be chosen. With local functions, the alteration 
of a control point results in an alteration of the 
constructed curve that is restricted to a local region 
about the point. The remaining regions are unaltered. 
Therefore local sections can be manipulated in an inde- 
pendent manner. The simplest local interpolants are 
the piecewise linear functions that do not vanish over, 
at most, two intervals defined by rj < rz < ... < r^-i- 
To further simplify the form of the multi surface trans- 
formation, the height of each interpolant shall be 
adjusted so that each interpolant integrates to unity. 
Then, the evaluation of the transformation at r^-i 
will reduce to by means of a telescopic collapse 

of terms. An explicit form of the normalized interpola- 
tion functions (Eiseman and Smith, 1980) is given for 
a case with the uniform partition (r^ = k, k = 1 , 2 , ..., 
N-l) 

1 2(2 - r) for 1 < r < 2 

( 11 ) 

0 for 2 < r < N - 1 


V r) 


0 

for 

1 < r < 

(r 

- k) + 

1 for 

(k 

- r) + 

1 for 

0 

for 

k + 1 < 


k - 1 

k - 1 < r < k 
k < r < k + 1 
r < N - 1 


( 12 ) 


1 0 for 1 < r < N - 2 

(13) 

2(r - N + 2 ) for N - 2 < r < N - 1 

Similarly, for the index value of i, a curve can 
be constructed 


F i <t> - c i, * L v^h.M -%} 

3=1 

1-1,2, N (14) 


where t is the curve parameterization and H^(t) is 
the integration of the multisurface interpolants. The 
tensor product form depends only upon Cjj and is 
given by J 

M-l 

T(r.t> - E,<r) + £2 H p (t) [e^ (r) - E ft (r)] (15) 

or alternatively by 

N-l 

T(r,t) - F,<t> + 22 G a (t) [ F a+l (t) “ F a (t) ] <16) 
a=l 

These two expressions are equivalent. The tensor 
product matches E-i or Fi at the extremities of 1 
and j. J 

When boundaries are to be specified, the corre- 
sponding data appear at the extremities of the values 
for r and t. Since the coordinate transformations 
are generally expressed in the form of a vector P(r,t) 
for the desired positions of all points in physical 
space, it is convenient to express the boundary specifi- 
cations in terms of the position vector. Thus, the 
boundaries are denoted by P(l,t), P( N-l , t ) , P(r,1), 
and P(r,M-l). To include the boundaries, the multisur- 
face transformation is performed again as above, but now 
with the actual boundaries Inserted. This results in a 
modification of T for both the r and t • directions. 
In each such directional construction, the actual bound- 
aries appear as end conditions for the corresponding 
variable whereas the remaining boundaries are solely 
generated by the control points. Thus, by subtracting 
T from the sum of both directional constructions . the 
actual boundaries become end conditions for each varia- 
ble. This process follows a Boolean sum format and 
upon simplification becomes 

Q(r.t) = T(r,t) + a,[l - G] (r)J ^P( I ,t) - F ] ( t>] 

+ a 2 G N _,(r)[p(N-l,t) - F N (t)] 

+ a 3 [l - H,(t>] [p(r,l) - E 1 (r)j 


+ a 4 H M 1 (t)[p(r,M-l) - E M (r)] 


(17) 
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where each of the four terms following the tensor prod- 
uct T(r,t) represents a transfinite conformity to a 
boundary when each ai Is 1. By setting any <x\ to 0, 
the corresponding boundary becomes available for 
free-form modeling by means of the control points. 

In the order listed, the boundaries are for r = 1, 
r » N - 1 , t » 1 , and t = M - 1. Further details are 
discussed by Eiseman (1987). 

RESULTS 

The basic interactive process of TURBO is illus- 
trated with a simple geometry in Fig. 2. It starts 
with construction of a simple control net and surface 
grid. An initial grid is generated and examined. If 
desired, the grid structure can be improved through the 
dynamic movement of the control points. In this example 
a control point was moved to obtain dense meshes in the 
middle of the flow region. 

Figure 3 shows an H-grid generated for a cambered 
turbine stator vane. Only those basic interactive fea- 
tures illustrated in Fig. 2 were used in this example. 
Figure 3(a) shows an initial control net; Fig. 3(b) 
shows an initial grid. By scrolling graphic images for 
one view (combination view of the control net and grid) 
to another (control net only or grid only) and by 
expanding the graphic image around a center of expan- 
sion, a user can closely examine the grid and determine 
what modification of the control net is needed to 
improve the grid structure. The modified control net 
and grid are shown in Fig. 3 parts (c) and (d), respec- 
tively. The modified grid has a densely clustered mesh 
around the leading edge and a more nearly orthogonal 
grid near the vane surfaces. 


boundary. Detailed stator vane geometry, cascade flow 
conditions, laser anemometer survey measurements, and 
vane surface static pressure measurements are presented 
by Goldman and Seasholz (1982). To operate the cascade 
facility, ambient air from the test cell was drawn 
through the cascade and exhausted into the laboratory 
altitude exhaust system. The test conditions in the 
cascade were set by controlling the pressure ratio 
across the vane row with two throttle valves located in 
the exhaust system. A hub static tap located downstream 
from the test section, where the flow was assumed to be 
nearly circumferentially uniform, was used to set this 
pressure ratio. The hub-static- to i nl et-total -pressure 
ratio was maintained at a value of 0.65. This corre- 
sponds to a static- to inlet-total-pressure ratio of 
0.687 at a 50-percent span where a two-dimensional 
inviscid computation was made. 

Core turbine stator vane geometry at the mean sec- 
tion is shown in Fig. 8(a). Also shown are three axial 
locations where measured data are compared with computed 
results. In Fig. 8(b), laser measurements are compared 
with computation at 10.2, 50, and 90.1 percent of the 
axial chord at the 50-percent span. The computed criti- 
cal velocity ratios (V/V cr ) agree well with the laser 
measurements. Figure 9 shows a good agreement between 
the blade surface measurement and computation. For this 
comparison the vane surface static pressure measurements 
are used to obtain the critical velocity ratios by using 
the following relation: 



Many grid generation programs do not allow local 
mesh control as TURBO does. For the TURBO program the 
initial grid may be generated by any program familiar 
to a user; then an initial control net may be obtained 
from the grid by attachment to essentially reproduce 
the grid structure. The initial control net shown in 
Fig. 4(a) was created by attachment. Figure 4(b) shows 
an initial grid of a compressor rotor blade. Once the 
control net is created, the interactive process to be 
followed is the same as described previously; the 
initial control net is modified to the one shown in 
Fig. 4(c), and a new grid (Fig. 4(d)) is generated from 
the modified control. A closer view of the initial and 
improved grid is presented in Fig. 5. The improved 
grid is more othogonal and has slope continuity across 
the periodic boundary. The shape of the control net 
can easily be changed by using an interactive process 
illustrated in Fig. 6. In part (a) of the figure, a 
user first chooses a control line to be changed and 
then picks a hinge point. In part (b), control point 1 
is moved to point 2 by moving the mouse of the work- 
station. Point 1' automatically moves to point 2' 
while maintaining a pitch with the points 1 and 2, 
respect! vely. In part (c), the desired shape of the 
control line is obtained by choosing a menu called 
"rubberbanding. " 

The finite-volume LU implicit scheme was used to 
simulate an inviscid flow in a turbine cascade on the 
grid shown in Fig. 7. Part (a) of the figure shows the 
control net from which the grid shown in part (b) was 
generated. Part (c) is a closer view of the grid. 

This grid is nearly orthogonal in the vicinity of the 
solid vane surface and has slope continuity across the 
periodic boundary. For this cambered vane, the skewness 
of this single block grid was reduced to some degree by 
using the "rubberbandi ng" feature and free-form periodic 


where p is static pressure and pt is total pressure. 
A static pressure contour plot is presented in Fig. 10 

without comparison. The pressure was normalized by 
using the inlet static pressure. 

CONCLUDING REMARKS 

The LU approximate factorization scheme was demon- 
strated to be robust and efficient in the transonic- 
flow and very high Mach-number-f low regimes. In these 
previous works, the main interest was in improving the 
efficiency of the scheme and therefore simple isolated 
airfoil and inlet geometries were used. In the work 
presented, the finite-volume LU implicit scheme was 
applied to obtain the solution of the inviscid flow in 
a turbine cascade on a grid generated by TURBO. This 
is a first step towards simulation of more complex 
flows in turbomachinery. 

TURBO is a menu-driven interactive algebraic 
grid generation program. It is being built around the 
control net which guides the coordinate curves. From 
the control points sparsely distributed over the flow 
domain, TURBO can easily modify the grid structure of a 
local region (a domain of dependence) without affecting 
the grid outside of it. For the complex flows in turbo- 
machinery, the control point formulation is a very use- 
ful method because of this precise local mesh control 
capability. With the advent of faster graphics worksta- 
tions with more memory and better linking to mainframe 
supercomputers , interactive grid generation by means of 
the control point formulation will become a more effec- 
tive and useful approach. 
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FIG. 1. - CONTROL NET - AN EXAMPLE IN TWO 
DIMENSIONS. 
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CONSTRUCT SIMPLE CONTROL NET; GENERATE 
SURFACE GRID. 




TRANSLATE CONTROL POINT (E.G., TO GET SELECT CONTROL POINT TO MODIFY GRID. 

FINE MESHES IN THE MIDDLE). 



RE-COMPUTE. EXAMINE NEW GRID - ACCEPT OR REPEAT 

PROCESS. 


FIG. 2. - BASIC INTERACTIVE PROCESS. 
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(a) INITIAL CONTROL NET. (b) INITIAL GRID. 

(C) MODIFIED CONTROL NET. (d) MODIFIED GRID. 

FIG. 3. - GRID GENERATION FOR A CAMBERED STATOR VANE. 





(b) INITIAL GRID. 




(C) MODIFIED CONTROL NET. 



(d) MODIFIED GRID. 

FIG. R. - GRID GENERATION FOR A COMPRESSOR BLADE. 
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(b) MODIFIED GRID. 


FIG. 5. - CLOSER VIEW OF THE INITIAL AND MODIFIED GRIDS. 
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(b) TRANSLATE A CONTROL POINT FROM 1 TO 2. 



(C) USE "RUBBERBANDING" FEATURE. 

FIG. 6. - CONTROL LINE MODIFICATION BY MEANS OF AN 
INTERACTIVE FEATURE. 
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(3) CONTROL NET. 



FIG. 7. - GRID GENERATION FOR A TURBINE 
CASCADE SIMULATION. 
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FLOW 


-*-y I 


SURVEY LOCATION AS 
PERCENT OF AXIAL CHORD 



(a) CORE TURBINE STATOR VANE GEOMETRY, SURVEY LOCATIONS, CIR- 
CUMFERENTIAL POSITION (y). 




FIG. 8. - COMPARISON OF COMPUTED RESULTS WITH EXPERI- 
MENTAL DATA. 
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NORMALIZED AXIAL DISTANCE, Z-Z LE /CHORD 
FIG. 9. - LOADING DIAGRAM. 
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INLET MACH 0.210 ALPHA 0.00 

PS/POIN ABS MIN 0.704 MAX 1.189 INC 

FIG. 10. - STATIC PRESSURE CONTOUR. 


0.034 


15 


2. Government Accession No. 


Report Documentation Page 

Space Administration 

1. Report No. 2. Government Accession No. 

NASA TM-101412 

4. Title and Subtitle 

Application of a Lower-Upper Implicit Scheme and an Interactive Grid 
Generation for Turbomachinery Flow Field Simulations 


3. Recipient’s Catalog No. 


5. Report Date 


6. Performing Organization Code 


7. Author(s) 

Yung K. Choo, Woo-Yung Soh, and Seokkwan Yoon 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


8. Performing Organization Report No. 

E-4374 

10. Work Unit No. 

505-62-21 

11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


| 15. Supplementary Notes 

Prepared for the 34th International Gas Turbine and Aeroengine Congress and Exposition, sponsored by the 
American Society of Mechanical Engineers, Toronto, Canada, June 4-8, 1989. Yung K. Choo, NASA Lewis 
Research Center; Woo-Yung Soh, Sverdrup Technology, Inc., NASA Lewis Research Center Group, Cleveland, 
Ohio 44135; Seokkwan Yoon, MCAT Institute, NASA Ames Research Center, Moffett Field, California 94035. 

16. Abstract 

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a turbine cascade. This 
approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, 
which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to 
reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to 
generate grids. This program uses the control point form of algebraic grid generation which uses a sparse 
collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct 
advantage of TURBO compared with other grid generation programs is that it allows the easy change of local 
mesh structure without affecting the grid outside the domain of dependence. Sample grids are generated by 
TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the 
LU implicit scheme on the grid generated by TURBO. 


17. Key Words (Suggested by Author(s)) 

Numerical scheme 
Grid generation 
Turbomachinery 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 02 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21 . No of pages 
16 


22. Price* 


NASA FORM 1626 OCT 86 


For sale by the National Technical Information Service, Springfield, Virginia 22161 





